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In this paper I discuss a mechanism for ferrimagnetism in (l+l)-dimensions. The mechanism is 
' related to a special class of interactions described by operators with non-zero Lorentz spin. Such 

operators are present in such problems as the problem of tunneling between Luttinger liquids and 
the problem of frustrated spin ladder. Exact solutions are presented for a representative class of 
models possessing a continuous isotopic symmetry. It is shown that the interactions (i) dynamically 
generate static oscillations with the wave vector dependent on the coupling constant, (ii) give rise 
to spontaneous breaking of this symmetry at T = accompanied by generation of the magnetic 
moment and appearence of gap less modes with a non-relativistic {ferromagnetic) dispersion E ^ k'^, 
(iii) generate massive (roton) modes. 
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, I. INTRODUCTION 
^ ■ 

. Thotigh general stability of critical points is determined by scaling dimensions of the perturbing operators, to 

i' ' determine the ultimate destination of the renormalization group (RG) flow is much more difficult problem. Availability 

''O ■ of non-perturbative methods makes this task easier in (1+1)- and two dimensions. 



Ch ' In (H-1)- or two dimensions critical points possess conformal symmetry with the Lorentz (0(2)) symmetry being 
its part. F] If the relevant perturbation preserves the Lorentz symmetry, the most it can do is to open a gap in a 



o 
o 



> 

00 



part of the spectrum. That is exactly what happens in all known solvable examples (for instance, in the sine-Gordon 
model) . 

It may happen however, that the perturbation violates the Lorentz symmetry. Such perturbations certainly make 
sense in condensed matter physics where preservation of the Lorentz symmetry is not required. There are quite 
' a few models of condensed matter physics which can be treated as critical models perturbed by relevant Lorentz- 
. symmetry-breaking operators. These operators have non-zero Lorentz spin. Since in (1-1-1)- or two dimensions the 
' Lorentz symmetry at criticality is extended to conformal symmetry, Lorentz spin is also called conformal spin. At 
criticality where all two-point correlation functions follow power law behaviour, an operator 0(z, z) is characterized 
by its scaling dimension d and conformal (or Lorentz) spin S defined as 

■ < 0(zi, zi)0(z2, -Z2) >= Az-^^+'^-z"^"-'^ (1) 

S ■ 

I ■ where z, z are holomorphic coordinates defined as z ~ t + ix, z — t ~ ix (in (l-l-l)-dimensions r is Matsubara time). 
■ O \ Since only bosonic operators can appear as perturbations, possible Lorentz spins are integer. Furthermore, since in 
^ unitary theories conformal dimensions are positive (d ± S* > 0) and d < 2 for relevant operators, this leaves us with 
only one choice: S* = ±1. 

The simplest example of = ±1 operator is a conserved charge. At criticality such perturbation generates incom- 
mensurability. As an illustration one can consider a change of the chemical potential for massless Dirac fermions: 
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A = J dTdx[R+dR + L+dL- ^i{R+R + L+L)] (2) 

This perturbation is removed by a simple transformation of the fields 

R^e^f'-'R, i->e-i^^i (3) 



^Stricktly speaking, the Lorentz symmetry is violated in the models with complex isotopic symmetry, where different sectors 
of the spectrum may have different velocities. However, this does not influence the argument made in the text. 
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and thus leads to a shift in the corresponding Fermi momentum and generates static oscillations with the new wave 
vector whose magnitude depends on the value of /z. The perturbed theory remains critical in the infrared (IR). 

The situation becomes less clear when the perturbation is not a constant of motion. Such non-trivial perturbation 
corresponds to a non-holomorphic operator, that is, in the notations of Eq.(0), an operator with d ^ ±S. 

A well studied example of 5 = ±1 perturbation is the the so-called Z^r 'chiral clock' model of statistical physics 
(see for review). In fact, this model is not too good for statistical physics purposes since its Boltzmann weights 
are not positively defined. However, the (l-l-l)-dimensional version of this model makes perfect sense and, as was 
demonstrated in Q, describes a tunneling between two spinless Luttinger liquids. It was established that for the 
chiral clock model the relevant S = ±1 perturbations generate incommensurability. 

An immediate generalization of the problem of tunneling includes spin. The tunneling between two spin-1/2 
Luttinger liquids is described by the Hamiltonian 

H = Hi'l + H^^l + I dx{T+ + T_), (4) 

T+ =t(i?+^i?2,. + i7.c), r^=t{L+^L2^, + H.c) (5) 

where i?i,2 and L12 are the right- and the left-moving fcrmions on chains 1 and 2; H^l describe interacting electrons on 
the individual chains. Due to the intra-chain interaction the fermionic operators acquire non-trivial scaling dimensions 
such that 

<T,(r,.)T,(0,0)>^^^^^^ (6) 

where 6 depends on the interaction (here I do not discuss the effects related to difference between the charge Vc and 
the spin Vg velocities of the excitations and set = = !)• According to definition operators T± carry conformal 
(Lorentz) spin 5 = ±1. 

Another problem much discussed in the literature is the problem of frustrated spin-1/2 two-leg Hcisenberg ladder 
(alias the zig-zag ladder). In the decoupled limit, two S=l/2 chains represent an SU {2)i x SU {2)i WZNW theory. 
Each SU{2)i WZNW model has its matrix field gi{x) {i = 1,2) = ei{x) + in^ • ct, (ni^2 are staggered magnetizations 
of chains 1 and 2 and £1,2(2;) are the staggered energy density operators). For the frustrated ladder the interchain 
interaction contains is dominated by the so-called twist term The most general SU(2)-invariant form of this term 
is 

Ctwist = Am ■ d^n2 + Beid^e2 + [1^2] (7) 

In the Heisenberg zigzag ladder the bare value of B is zero but it is generated in the course of RG (due to the fusion of 
the A-operator with scalar marginal ones), where the leading interchain interaction gives the following contribution 
to the action density j^: For a single chain the two-point function of the operators n, e decays as 



which means that 



< n(T, x)n(0, 0) >=< e(T, a;)e(0, 0) (a;^ + t^)-!/^ (8) 



< 0(t, x)O(0, 0) {x" + t2)-i [(t - i.x)-2 + (r + \x)-^] (9) 



According to definition (^, this operator is a sum of two operators with d = 2,5 = ±1. 

Using the procedure suggested in ||^, one can reformulate the zig-zag ladder model as a model of four Majorana 
fermions with the Lagrangian density given by 

^ = -^0 + C-twist, 

3 1 I 
^0 = ^[^'X.aidr - '^^^x)Xa + -jXa{dr + ''^^x)Xa] 

Ctwist = {51X1X2X3X0 + 52[xoXiX2X3 + X3X0X1X2 + X2X3X0X1] + [(X ^ X)]} (10) 

where (71,(72 are linear combinations of A,B. In the limit A — B we get the perfect 0(4) symmetry gi — §2, with a 
suitable choice of A, B one can get §2 — which corresponds to the model I consider later in the paper. 

Operators with non-zero conformal spin may also appear in models describing tunneling between edges of incom- 
pressible Fractional Quantum Hall (FQH) states 0. 
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Since operators with non-zero conformal spin break Lorentz symmetry, it lifts the restrictions on the low energy 
spectrum. With the Lorentz symmetry being preserved one has two choices: either to have a linear spectrum or 
to have a spectral gap. Without Lorentz symmetry one can imagine a quadratic gapless spectrum which would be 
naturally associated with a ferromagnetic Goldstone mode. Thus, if a model in question has a continuous isotopic 
symmetry, the presence of 5 = ±1 perturbation may lead to a spontaneously breaking of this symmetry at T = 
with a formation of a net 'magnetic' moment and the above-mentioned quadratic spectrum. 

In all examples discussed above (except the one involving the edge states) the perturbations with S* = ±1 appear 
together. This introduces an additional complification because a fusion of two operators with opposite Lorentz spins 
generates under RG a relevant operator with zero spin. However, in a certain range of parameters of model this 
secondary flow does not catch up with the flow of the original coupling constant (for model (^ this is always the case) 
0' [H' H] (^^6 also [Q) where this topic is extensively discussed). This gives room for an energy scale on which the 
effect of interaction between the operators with S — 1 and S — —1 can be neglected. On this energy scale one can 
drop one operator and consider a perturbation with one sign of conformal spin, for example S' = 1 |p^ . In FQH 
systems, where excitations are chiral, appearence of a single 5 = ±1 perturbation is more likely. 



II. A SOLVABLE MODEL 



In the subsequent sections I discuss a series of exactly solvable models which realize the symmetry breaking scenario 
outlined above. These models are related to Wess-Zumino-Novikov-Witten (WZNW) model. The latter model 
together with the related lattice models has been used extensively to address the problems in magnetism (see, for 
example, [0, [Q, (l^). Some of these lattice models are rather similar to the zig-zag spin ladder (p^. 

The particular model I discuss is described by the following action: 

A^W^ + haJd''x:.r{x)^aAx)- (11) 

where is the action of the critical Wess-Zumino-Novikov-Witten (WZNW) model on the SU(2) group with level 
k > 2, $a,a are the primary field from the adjoint representation, Ja are the left Kac-Moody currents and ha is a 
constant vector. The restriction on the value of k is related to the fact that for fc = 1 field ^aa does not exist. The 
suggested solution can be easily generalized for other symmetry groups and coset models. 

To make the discussion self-contained, I recall several basic facts about WZNW models. A WZNW model describes 
a matrix field g{x) defined on a group G whose dynamics is governed by the action 

Wig) = Y^J d^xTvid.g-'d^g) + kr[g] (12) 

r[.9] = J d^xe'^f'^Trig-'d^gg-'di.gg-'d^g) 

The action contains two parameters - the coupling constant A and integer number k. The model has a global Gr x Gl 
symmetry being invariant under the transformations 

g{x) ^ Vg{x)U (13) 

Hence there are two conserved charges associated with left and right global shifts of the matrix field g. In the perturbed 
model ( |rT| ) the right symmetry is broken, but the left symmetry is not. As we shall see below, this symmetry is broken 
spontaneously at zero temperature. 

Model ( |l^ ) has a stable critical point; the critical value of the coupling constant is A* = k^^. In what follows I 
shall distinguish between WZNW model as such and critical WZNW model described by action ( [l^ ) with A = k~^. 

At the critical point WZNW models possesses a higher symmetry: its right (left) currents become holomorphic 
(antiholomorphic) : 

J=^g-^dg, 9J = 0, 

J = -^5%^\ dJ^O (14) 
(where d = dz, B ^ dg] z — t + ix, z = t — ix) and satisfy the Kac-Moody algebra: 
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[J-{x), J\y)] = '^5'^'5'{x -y) + ir'^r{y)5{x - y) 



(15) 



The critical WZNW model can be conveniently described using the Hamiltonian formalism: 

H = J2 / J°(a;)J"(a;) : + : J'^(x)J"(a;) :] 

^ a— 1 

where Cy is the quadratic Casimir in the adjoint representation. For the SU(2) group = 2. 

At criticality any field which is local in g can be represented as a linear combination of mutually local operators 
composing a basis in the operator space. This basis of fields contains the primary fields ^^^'^ and their descendants. 
The primary fields transform under Gr x Gl as tensors belonging to irreducible representations of the group; the 
descendants are generated from the primary fields by fusion with the Kac-Moody currents. Primary fields from higher 
representations can be generated by fusion of the fields from the fundamental representation. However, this process 
terminates after certain number of fusions depending on the value of k. 

To make the discussion more concrete, I concentrate on the SU(2) group. In this case the primary fields are 
characterized by the value of spin j — 1/2, 1, 3/2, .... The field with the smallest spin j = 1/2 is the g-matrix itself. A 
fusion of two g-operators generates a field with j — I- For the SU(2) group this representation is isomorphic to the 
adjoint one. The maximal spin one can achieve by fusing (7-matrices is equal to fc/2. Therefore for k = I model there 
is no primary field with j — 1 and model (pi] ) is not well defined. 



III. BETHE ANSATZ SOLUTION 

To solve model ( p] ) I use the Bethe ansatz solution of WZNW model |l6). This solution describes not just the 
critical point, but the entire RG flow towards it. As was shown in ||l^, in the vicinity of the critical point the leading 
irrelevant operator responsible for this flow is 

V = X:JaJa{x)<i>a,a{x): (16) 



In the presence of the right magnetic field Hr this operator generates perturbation (jll]). Indeed, this field is coupled 
only to the currents of right chirality Jg. In a finite magnetic field these currents acquire a finite expectation value 

< J;, >=(5'^''(ij(«)/27r) (17) 

and therefore in the leading order in H the irrelevant operator (^6|) is transformed into a relevant perturbation ([ll| ) 
with 

^ XiHi'^^ /2n) (18) 

Thermodynamic Bethe Ansatz (TBA) equations for SU(2) WZNW model of level k in magnetic field have the 
following form |l^ ]: 

„(v) = Ts * ln[l + e<^"+i('')/'r][l + e^.-iC'')/^] _ mS^^e^''/'^ - m5n,ke"''"''^ , (19) 



e 



lim f!i = -Hr, hm —^Hl 

n — ^ — oo 71 n — 'oo Tl 

where n = — oo, ...oo. Hr and are right and left 'magnetic' fields corresponding to the two conserved charges. In 
accordance with the above discussion, I will keep only the leading contribution in Hr = H and keep infinidecimal. 
The free energy is given by 

F/L ^ ~mT J dvie''''/^Sn,o + 5n,ke~''''''^) ln[l + e^-t")/"^] (20) 

[L is the system's size) and 

oc 



S * f{v) = / du- ^^""^ 



4cosh[7r(?; — u)/2] 
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We first study these equations at T = 0. In tiiis case all e„ with negative n are of order of nH, €q{v) is positive at 
V — oo and changes its sign at B — \a{m/ H), e„ with positive n are also positive except of ek{v) which changes 
its sign from negative to positive at some value Q > 0. The ground state energy can be expressed in terms of eo and 
Efc which satisfy the following integral equations: 



AuD{v - u)ek(u) = -me-"""'^ 



duK{v ~ u)eo(u) 



duD{v - u)eo(u) = -rne'^''/^ + H + duK{v - u)ek{u) 



where the Fourier transforms of the kernels are 

„, , tanhwe'^l'^l 

D{uj) = —-, — 

2 sinh ku 



tanho; 
2 sinh kuj 



(21) 
(22) 

(23) 



The limits B and Q are determined by the conditions eo{—B) = 0, ek{Q) = 0. 

To determine the spectrum, we also need to know the distribution densities po, Pk which satisfy similar equations: 



^ /2 



duD{v — u)pk{u) ~ — e 



duD{v — u)po{u) = — e 



oo 
oo 



(24) 
(25) 

(26) 



Since we are interested only in the leading asymptotics in H, we can neglect in Eq. (|22|) (and respectively pk in 
Eq.(p5[)) and solve these equations as Wiener-Hopf ones. Eqs.(pl|, p4| ) in this case also become Wiener- Hopf equations. 
The case fc = 2 is somewhat special and will be treated separately. For fc > 2 the solutions at large Q, B (small H/m) 
are given by 



duK{v — u)po{u) 
duK{v — u)pk{u) 

These distribution densities determine the conserved charges (right and left 'magnetic' moments): 

1 [°° 1 [^^ 

Ql = t; I <^upk{u), Qr^ - dupo{u) 

^ Jo ^ J-oo 



eo(w) 
efc(w) 



-r^ 27r 



■kH 



duj 
2^ 



-\ui{v-Q) (+) , 



e -6r^(c.), 6^(0,) 



icj(-7r + 2icj)G(+) {uj)Gi+){Q) ' 
7r(fc - 2)me^'^'3/2 
' (tt - 2iu;)(7r - i/cw)G(+) (w)G(-) (-i7r/2) ' 



pW(c.) = 4+)(^) ("~^^"; 



where 



TOexp(-7rB/2) = H 



G(+)(i7r/2) 



G(+)(0) ' 
exp(-^g/2) = PkiH/m)^, 
k tan(7r/fc) 



(27) 
(28) 
(29) 

(30) 
(31) 



27r(fc + 2)[G(+)(i7r/fc)]2 



|fe/(fe-2) 


■G(+)(i7r/2)" 




G(+)(0) 



and 



lujk + 



r(i + i|^)r(i + ^) 
V2^r(i + ^f^) 



(32) 



with G(-)(w) = GW(-w) and D{u) ^ G^-'> {uj)G^+'> (uj). 
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To obtain the modified TBA equations describing model (^Tj) I simply replace cq in Eq. (|T^) by its zero temperature 
approximate value (p7|). The result is 

en{v) = Ts * ln[l + e^"+i('')/^][l + e^"-i(")/^] + 6n,is * eo{v) - m,5„,fce--"/2, {n = 1, ...) (33) 
lim — = Hl 

n— >oo n 

The free energy is then given by 

F/L = -mT J dve'^"/^ ln[l + e<^'=(")/'^] (34) 



Let us return to T = solution. From equations ( |2^, 29 ) we can extract the following information. First of all, I 
observe that approaches zero at two points: at v = Q and v —> oo. The latter point I identify with zero momentum; 
then point v = Q corresponds to the wave vector 

f°° f ,1 (fc + 2) 

(5P = 27r/ dvpk{v) ^2t:p[^'{lu ^0) ^m'yk{H/m)T^ (35) 

where 

Ik 



G(+)(i7r/2)G(+)(0) 



This incommensurate wave vector scales with H exactly as one expects it to scale (that is 5P ~ iJ^/^^^'*)), taking 
into account that the scaling dimension of the perturbation in Eq. (pT|) is d = 1 + 4/(fc + 2). The excitation spectrum 
in the vicinity of w = Q is linear; the velocity is given by 

V 1 deu{v) (fc-2) 

^^=2vrp,(g) dv 1"-^ = ^^^ (3^) 

where Vl is the velocity of the left-moving particles (I have been working in the system of units where Vl = 1). 

The fact that approaches zero also at w ^ oo means that there are soft modes at zero momentum. Since according 
to Eqs.(28, 2^) efc('^), Pfc(w) behave as at small uj, their real space asymptotics at t; — > oo are £k ^ Pk ^ v~'^ . Since 
the momentum is given by 

/•oc 

P{v)=2tt dvpk{v) (37) 

J V 

this means that the spectrum is quadratic. Using this formula and Eq.(p9|) we find the dispersion law: 

The maximum of the energy is reached at P ~ 6P and is of order of 6P. 

For k = 2 the formulae for the velocity and the dispersion law should be modified: 

Q = 27r/g+iln(4g), g ^ [HG^+Hi7r/2)/m]^ (39) 

Vr/Vl^^9-H^, (40) 

'^^ ^ H[Gi^Hin/2W «-P{-[-™/^G(+)(iV2)]n (41) 



Taking T limit in Eqs.(|3^) I obtain the following expressions for the energies: 

er,(w) = 

sinhncj = Q (+) sinh(fc - n)^ (+) 

sinhfcc.'' ''^ ^""'^ sinhfco; ^ [ 
n < k 



(42) 



e„+fe(c.) = e-"l"lei+H..)ei"'3 (43) 
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Let us study asymptotics of these energy functions. In real space we have 



1 



du 



u) 



(44) 



From this expression it follows that £n+k{v) decay as u ^ on both infinities. On the other hand, en{v) with n < k 
decay as v^"^ only at v >> Q. At v < Q they have roton-likc minima. For Q — v >> 1 I obtain from Eq.(|4^): 



« k-^sm{Tm/k) \e+{m/k)e'''-^-Q^^^ + 4{m/k)c-<^+^^^^' 
^ 2fc-i sin(7rn/fc)[e+(i7r//c)e+(i7r/fc)]i/2e-'^(Q+^)/2fe cosh[7r(t; - VQ)/k] 

which corresponds to the relativistic spectrum with spectral gaps for k ^ 2 {n ~ 1, ...A; — 1) given by 

"2(fc-2)cot(7r/fc)1^/^ 



Mn = Mq sin(7rn/fc), Mq = 



nk{k + 2) 



6P 



For k ~ 2 the gap is exponential in i7 ^. The energy minimum for e„ with n = 1, ...fc — 1 occurs at 

_ k j (fc - 2)tan(7r/fc) 

" ^ " 2^ \47r(fc + 2)[G(+)(iVfc)]2 



(45) 



(46) 



(47) 



This means two things. The first one is that the massive modes are centered at incommensurate wave vector. The 
second is that in the rapidity space they are very close to Q. 



IV. BETHE ANSATZ DERIVATION OF THE LOW ENERGY EFFECTIVE ACTION 



Thus we have the following regions in rapidity space where low-energy excitations are located: 

(i) at V » Q there are gapless modes e„ (n = 1, ...) with the spectrum ~ 

(ii) at V << Q there are gapless modes e„ {n = k + 1, ...) with the spectrum ~ 

(iii) at V — Q there is a gapless mode Ck with the spectrum ~ {v — Q)', 

(iv) at V K Q there are massive modes e„ (n = 1, ...fc — 1). The massive modes can be treated as low energy 
excitations only for fc >> 1 when their masses are much smaller than SP. 

In the momentum space the spectrum in sectors (i),(ii) is quadratic e(P) ^ (see Figs. 1,2), the spectrum in 
sector (iii) is linear e(P) ~ |P— P{Q)\ and the spectrum in sector (iv) is massive relativistic (see Fig. 3). 
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^ p 

FIG. 2. Dispersion e„{P) (n = 1, ...k - 1) 




Let us first consider the case of moderate k when one does not need to consider roton modes. Then the low energy 
sector is described by the truely gapless modes (i - iii) . Since their positions in rapidity space are well separated from 
each other and the integration kernels in Eqs.(p3|) decay exponentially, one can derive separate sets of TBA equations 
for each mode. Comparing these equations with TBA equations for known integrable models one can deduce the 
effective action for the low energy sector. 

Let us derive TBA equations for sector (i). In this sector all energies are positive. Therefore it is convenient to 
have TBA equations in such a form where the integral kernels act on functions ln(l + e~'^"/^). Such equations would 
provide a ready T — limit. To get such form of TBA. The resulting equations read 

Tln[l + e^"(")/^]- A„„,*rin[l + e-^"(")/^] = A^^ifc*A„,fc*ei^°^(z;)+ni7L, n = l,..., v»Q (48) 

where €f\v) is the solution at T = and is given by Eq. (psf) . The kernels have the following standard Fourier 
transforms: 

A-nmi^) — coth |ti;| [exp(— |n — — exp(— |n + m||ti;|)] (49) 

These equations are valid for temperatures T << 6P - the maximum value of el^\v). In this region These equations 
coincide with the right chiral sector of the spin-(k/2) integrable ferromagnet. The latter ones can be extracted from 
TBA equations obtained in |jl^. In a similar fashion I obtain TBA for sector (ii): 

rin[l + e^"+''(")/^]-A„™*rin[l + e-^"+'=(")/^] =a„*e[^°^+niJL, n^l,..., v«Q (50) 

These equations describe the left chiral sector of the spin-1/2 integrable ferromagnet where the dispersion law is 
e{k) ~ k'^d{k). Therefore the parity in the ferromagnetic sector is broken. 

Another soft mode in the right chiral sector has the linear spectrum with velocity ( |3^ ) and is centered at the the 
momentum dP. Using the standard manipulations with TBA equations one can find its contribution to the specific 
heat and establish that this mode carries central charge 1. Therefore it is described by the non- chiral Gaussian model. 
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V. CONNECTION TO THE ZIG-ZAG CHAIN AT K = 2 



For k ~ 2 model ( [Til) simplifies considerably. In this case the critical WZNW action is equivalent to the theory of 
three massless Majorana fermions Xa, Xa (a = 1, 2, 3) and the adjoint operator and the currents are given by [|l5t 



^ab = XaXb, J" = XbXc 



Therefore for fc = 2 one can rewrite the perturbed action (|ll|) solely in terms of Majorana fermions: 
A-- 



d^x 



{Or + ldx)Xa + haXlX2X3Xa 



(51) 



(52) 



In this form the model closely resembles the model for the zig-zag ladder (|T^) (see also ||l^) and att 92 = it even 
coincides with model (|l^) exactly solvable by the Bethe ansatz. 
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